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4 Abstract. We perform a two-dimensional simulation by using an elec- 

5 tromagnetic hybrid code to study the formation of slow-mode shocks in col- 
e lisionless magnetic reconnection in low beta plasmas, and we focus on the 

? relation between the formation of slow shocks and the ion temperature anisotropy 

a enhanced at the shock downstream region. It is known that as magnetic re- 

9 connection develops, the parallel temperature along the magnetic field be- 

10 comes large in association with the anisotropic PSBL (plasma sheet bound- 

11 ary layer) ion beams, and this temperature anisotropy has a tendency to sup- 

12 press the formation of slow shocks. Based on our simulation result, we found 

13 that the slow shock formation is suppressed due to the large temperature anisotropy 
« near the X-type region, but the ion temperature anisotropy relaxes with in- 

15 creasing the distance from the magnetic neutral point. As a result, two pairs 

16 of current structures, which are the strong evidence of dissipation of mag- 

17 netic field in slow shocks, are formed at the distance |x| > 115 Ai from 
is the neutral point. 
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1. Introduction 

19 In the Earth's magnetotail, magnetic reconnection plays an important role in the con- 

20 version of magnetic field energy in two lobes into kinetic and thermal energy of plasmas in 

21 a plasma sheet. Since Petschek [1964] proposed the necessity of pairs of slow-mode shocks 

22 attached to the diffusion region in magnetic reconnection in order to achieve its efficient 

23 energy conversion rate, many studies have been devoted to this model. As for numerical 

24 studies, many MHD simulations have confirmed the existence of slow shocks along the 

25 reconnection layer. Namely, slow-mode waves, propagating from the neutral point toward 

26 two lobes and along the outflow jets at one time, steepen and result in the steady state 

27 two pairs of slow shocks as is suggested by Petscheck [Sato and Hayashi, 1979; Scholer, 

28 1989; Abe and Hoshino, 2001]. In addition, observations of ISEE [Feldman et al, 1985] 

29 and Geotail [Saito et al, 1995; Seon et al, 1995] satellites have shown the existence of 

30 such slow shocks in the Earth's magnetotail. Especially, Saito et al. [1995] showed that 

31 the variation in the ion temperature is much larger than the electron temperature varia- 

32 tion across slow shocks. This suggests that the ion scale dissipation mechanism would be 

33 strongly related to the formation of slow shocks in collisionless plasmas. 

34 Until now, however, there is no clear consensus on the formation of slow shocks in mag- 

35 netic reconnection by both hybrid and full-particle simulations, even though the formation 
so of slow shocks itself has been demonstrated by slow shock simulations without magnetic 
37 reconnection [Omidi and Winske, 1989; Karimabadi, 1995; Omidi et al, 1995], and by 
33 a Riemann problem of slow shocks [Fujimoto and Nakamura, 1994; Lin and Lee, 1995; 
39 Liu et al, 2011]. Attempts to investigate slow shocks in a large scale reconnection with 
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40 two-dimensional hybrid codes were first done by Krauss-Varban and Omidi [1995]. They 

41 showed a fine structure of the reconnection layer where the fast plasma flow is generated. 

42 They concluded that these transition layers did not confirm the properties of the expected 

43 slow shocks, and mentioned that the reason for this is due to the fact that the ion dis- 

44 sipation scale is comparable to the thickness of the developing plasma sheet. Lin and 

45 Swift [1996] and Lottermoser et al. [1998] also performed large scale hybrid simulations 

46 to investigate slow shocks in magnetic reconnection. Lin and Swift [1996] suggested that 

47 the isotropic Rankine-Hugoniot (RH) jump conditions of slow shocks were better satisfied 

48 with increasing distance from the neutral point in the case of no-guide field reconnection. 

49 Lottermoser et al. [1998] also showed the slow shock-like discontinuities, and suggested 
so that the downstream ions were not directly heated by slow shocks as was shown by Lin 
si and Swift [1996] but heated by the stochastic motion of ions. They discussed that the thin 

52 current sheet formed after reconnection became turbulent and such turbulent structures 

53 caused stochastic motion of ions. 

54 That being the case, what causes such a discrepancy of the formation of slow shocks 

55 between MHD [Sato and Hayashi, 1979] and kinetic treatments of plasmas? We suggest 

56 that an important kinetic modification in magnetic reconnection would be ion temperature 
s? anisotropy along reconnection layers. Observations by ISEE-3 and Geotail satellites in 
sa the Earth's magnetotail reported that the ion temperature parallel to the magnetic field 
so (Ti t \\) is higher than that perpendicular to the magnetic field (Tl,±) at the downstream 

60 region of slow shocks [Cowley et al, 1984; Hoshino et al, 2000]. Such ion temperature 

61 anisotropy is known to be produced by two plasma components, i.e., the convecting cold 

62 lobe ions with the velocity = cExB/B 2 and the PSBL (plasma sheet boundary layer) 
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63 ion beams which are accelerated and ejected from the diffusion region. Preceding studies 

64 for magnetic reconnection by kinetic simulations also support the fact that T^/T^ > 1 

65 at the downstream of discontinuities along the reconnection layer [Hoshino et al, 1998; 

66 Lottermoser et al, 1998]. On the other hand, theoretical studies about slow shocks with 
e? the one-dimensional RH relations suggest that conditions required to satisfy slow shock 

68 RH relations become restricted by the ion temperature anisotropy, T it \\/T^± > 1, in the 

69 downstream region [Lyu and Kan, 1986; Karimabadi et al, 1995]. They showed that slow 

70 shock solutions in anisotropic plasmas are greatly affected by the downstream tempera- 

71 ture anisotropy. If T- h \\/T- h ±_ > 1 in the downstream region, slow shock solutions can exist 

72 only for the limited upstream Mach number regime. This means that the enhancement of 

73 ion temperature anisotropy along the reconnection layer makes it harder for slow shocks 

74 to exist in magnetic reconnection. In this study, we focus on this point and perform 

75 kinetic simulations for magnetic reconnection. Then, the relation between such ion tern- 
re perature anisotropy and formation of slow shocks in magnetic reconnection is discussed 

77 in accordance with anisotropic RH relations. 

78 In the following sections, we will first refer to the simulation model; second, show 

79 results of the kinetic simulation for magnetic reconnection, and lastly discuss the nature 
so of discontinuities formed along reconnection layers in detail. 

2. Simulation Model 

si In our study, the two-dimensional electromagnetic hybrid code, in which ions are treated 

82 as particles while electrons as a mass-less fluid, is used to investigate the formation of slow 

83 shocks in magnetic reconnection. Algorithm of our hybrid code is based on the method 

84 of general predictor-corrector loops [Harned, 1982; Winske and Leroy, 1984]. Within our 
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as simulation, basic equations are as follows: 

— = -cV x E, (2) 

E = -Vp c - -V e x B + rjJ, (3) 

qirii c 

V ' = V -4^ VXB ' < 4 > 



Vi/i(vi)tZvi 

Vi = J -f§o , (5) 

/i(vi)dvi 



^ + (V e • V) Pe + 1Pc (V • V e ) - ( 7 - 1) vJ 2 = 0, (6) 

86 where is the velocity distribution function of ions and a charge neutral condition, 

s? qifii — en c = 0, is assumed. 7 is the adiabatic index and set to be 5/3 in this paper, 

ss As for electrons, we assume that the electron gas is isotropic and both the electron heat 

as flux and the viscosity stress tensor are neglected. In our simulation, various parameters 

oo are normalized by the parameters in the initial lobe (e.g., n = n/n , B = B/£> , and 

91 V = V /Vao, where no, B , and Vao are respectively the density, the magnetic field, and 

92 the Alfven velocity in the initial lobe). 

In addition, the spatial profile of the resistivity is given by 



7} (x, y) = i] + i] c cosh 



-2 



x\ 2 ( y 



2' 



(7) 



93 where i]q is the background resistivity due to ion-electron interactions. Here, the back- 

94 ground resistive length A r0 = rj c 2 / (AtiVao) is se t to be 10~ 4 A ; , where A; is the ion inertial 

95 length in the initial lobe. r\ c is the anomalous resistivity due to some instabilities in the 

96 diffusion region, e.g., the lower hybrid drift and the drift kink instabilities. This resistive 

97 term is assumed to be independent of time and its resistive length is equal to 3.6 x Ai. 
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98 l x and ly are the characteristic lengths to determine the size of the anomalous (electron) 

99 diffusion region and are set to be l x — 1.0 A ; and l y = 0.5 A i? respectively. 

The whole size of the two-dimensional system is —342 \ x < L x < 342 A; and —32 A; < 
L y < 96 A;. Grid intervals, A x and A y , are both equal to 1/3 A ; . Initially, a double Harris 
equilibrium is assumed and periodic boundary conditions in both x- and y-direction are 
imposed. The initial spatial profile of magnetic field is given as 



where y' c = 64.0 A; and the half thickness of an initial current sheet, 5 y , is set to be 1.2 Ai. 
The density is given by 



100 The ratio of the density at the center of the initial current sheet to the background one, 

101 i.e., n c /n , is set to be 4. We assume the uniform electron temperature, T e = 5 x 10~ 3 

102 (/? ej o = 10~ 2 in the initial lobes). Ions consist of two components: The current sheet ions 

103 and the background ions. The current sheet ion temperature is set to be T i c = 0.12, which 

104 satisfies the relation, n c (T itC + T e ) = Bq/(8tt). As for the background ion temperature, 

105 we assume T ij0 = 5 x 10 -4 (/3j i0 = 10~ 3 in the initial lobes). Both the current sheet and 
we the background ion temperatures are given isotropically. We initially input 160 super- 
107 particles per cell at the center of the current sheet (i.e., 32 super-particles per cell at two 



3. Results 

109 First we shall show the whole structure of reconnection, focusing on discontinuities 
no formed along reconnection layers. Next, the nature of such discontinuities in kinetic plas- 







108 



lobes). 
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in mas are investigated in accordance with MHD Rankine-Hugoniot relations in anisotropic 

112 plasmas. 

U3 Structure of reconnection layer: 

n4 In Figure 1, we show the whole structure of reconnection at time t = 415 fi; -1 . Magnetic 

us field lines continuously reconnect with each other at the center of the simulation box and 

He two pairs of reconnection layers are formed. In the region \x\ < 70 A;, ions experience 

n7 Speiser-type trajectories [Speiser, 1965; Nakamura et ai, 1998] and form a thin current 

us sheet whose half thickness is about an ion inertial length. In 70 A ; < \x\ < 110 A i; the 

n9 current structure becomes in some degree turbulent and its half thickness reaches about 

120 5 Ai. Then, the current appears to be concentrated in the areas along the PSBL as is 

121 predicted by MHD simulations. 

122 Figure 2 shows enlarged views of the reconnection layers at x < 0. From the top to 

123 the bottom, magnetic field lines, the out of plane magnetic field B z , flow vectors, the 

12 4 mean ion temperature (7]) = (Ty + 2T i)J _)/3, and the ion temperature ratio T it \\/T i; ± are 

125 shown. Vertical red dashed lines separate two regions: Region 1 and Region 2 as indicated 

126 in Figure 2(b). In Region 1, the reconnection layer reaches almost steady state and its 

127 global structure does not change over time except for small scale turbulent structures. In 

128 Region 2, the reconnection jet encounters the preceding plasmoid and the plasma flows 

129 diverge. At this time, the region \x\ < 150 A; is filled with plasmas which originate from 

130 two lobes. 

131 As is seen in Figure 2(a), magnetic field lines begin to bend from x ~ — 115 A; and pile up 

132 from x ~ —125 A ; . The out of plane magnetic field, B z , shows a clear quadrupole signature 

133 at —50 A; < x [Hesse and Winske, 1994; Nakamura et ai, 1998]. In —120 A; < x < —50 Ai, 
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134 the oscillations of B z appear near the central plasma sheet (CPS) [Karimabadi et al, 1999]. 

13 5 The outflow velocity at y = is ~ 0.7 Vao in Region 1 and ~ 0.5 Vao in Region 2. In 
us the piled-up region (x < —125 A ; ), ions are heated and the ion temperature anisotropy 
137 considerably decreases as is seen in Figure 2(d) and 2(e), while in —120 Ai < x (Region 
us 1) the ion temperature anisotropy remains high especially in the transition region (the 

139 region between upstream and downstream). The ratio Ti t \\/T^±_ is about 4-6 and in good 

140 agreement with Lin and Swift [1996] and Lottermoser et al. [1998]. 

141 Figures 3(a)-(c) show three different ion velocity distribution functions f(v x ,v y ) in the 

142 upstream, transition, and downstream regions, respectively. These ion velocity distribu- 

143 tion functions are constructed by using super-particles within the white squares shown in 

144 Figure 2(e). In the upstream region, only cold lobe ions exist with the E x B drift velocity 

145 V y ~ cE z B x /B 2 ~ —0.15 Vao- In the transition region, both cold lobe ions and PSBL 
«6 beam ions exist. Characteristic speed of the PSBL beam ions is VJbeam ~ 1-2 Vao, an d it is 
147 known that these PSBL ions are accelerated in and around the diffusion region [Hoshino 
us et al, 1998]. Here, note that the ion temperature evaluated from these two components 
149 in Figure 3(b) gives > T\ t ± in the transition region. In the downstream region, ions 
iso are considerably heated, and the shifted-Maxwellian distribution with the outflow velocity 

151 V x ~ —0.7 Vao is observed. 

152 Next, to investigate fine structures along reconnection layers in more detail, Figures 4-6 

153 respectively show the cross-sectional views of discontinuities at x — —88.3 Ai, —115.0 A; 

154 and —145.0 A; of Figure 2. (These locations are indicated by black arrows at the bottom of 

155 Figure 2(e).) The horizontal axes correspond to y axes of Figure 2, which are normal to the 
ise initial current layer. From the left top to the right bottom, the ion density, bulk velocities 
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157 V x and Vy, the total magnetic field strength, x and z components of the magnetic field, the 

158 current along the initial current, the ion temperature ratio T^/T it ±, and the anisotropic 

159 parameter e = 1 — (f3\\ — (3±)/2 are shown, respectively. The anisotropic parameter e is 
mo useful to discuss the net effects of the temperature anisotropy. The vertical dash-dotted 
mi line stands for the boundary between the upstream and transition regions. The vertical 

162 dotted line stands for the boundary between the transition and downstream regions. The 

163 main judgmental standard points to determine these boundaries are the changes of both 

164 outflow and inflow bulk velocities. Note that the horizontal dashed lines shown in both 
us T i; \\/T it ±- and e-plots stand for the isotropic temperature baselines. 

166 At x — —88.3 A;, the current J z concentrates in the CPS, even though slow shock-like 

167 discontinuities are formed along the reconnection layers. At x = —115 A;, the enhancement 

168 of J z is seen around the transition region. Such bifurcated structures of the current 

169 indicate that most of the magnetic field energy is converted into kinetic and thermal 
™ energy of plasmas not in the CPS but near the transition regions. As for the magnetic 

171 field, the out of plane magnetic field B z is confined to the transition region. The changes 

172 of the density, the bulk velocity, and the magnetic field also show the slow shocks-like 

173 behavior across these discontinuities. At x — —145 A i5 a pair of current layers can be 

174 seen clearly. Note that the downstream of the pair of the current layers is the plasmoid, 

175 and that the inherent heating mechanism would be different from that in the quasi-steady 

176 state region (Region 1). In Region 2, since the radius of curvature of the magnetic field at 
I?? the CPS becomes as large as the ion gyro-radius, the ion motion becomes stochastic and 

178 would contribute to the plasma heating [Bilechner and Zelenyi, 1989; Lottermoser et al, 

179 1998]. 
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iso As is shown in Figures 4-6, Tj^ is larger than T- h ± in the transition region. The ion 

181 temperature anisotropies are reduced across the transition region, but we can observe a 

182 finite ion temperature anisotropy in the downstream region for x = —88.3 A;. In the 

183 downstream region, since plasma betas are larger than those in the transition region, 

184 net anisotropic effects are significant and may affect the structure of discontinuities and 

185 the plasma instability. Also it should be noted from e-plots in Figures 4-6 that this net 

186 anisotropic effect of ions becomes large with decreasing the distance from the neutral 
is? point. In —100 A; < x, the downstream temperature anisotropy is large and the fire-hose 

188 unstable condition is often satisfied as is shown in the e-plot of Figure 4. (Note that the 

189 fire-hose unstable condition is p\\ — p± > -B 2 /(47r), i.e., e < 0.) At x — —115 A; (Region 
wo 1) and at x — —145 A; (Region 2), the ion temperatures are almost isotropic. From 

191 the point of view of RH relations in anisotropic plasmas, such downstream temperature 

192 anisotropy is important to discuss the discontinuities [Lyu and Kan, 1986; Karimabadi 

193 et al, 1995]. This effect in magnetic reconnection is discussed in the following subsection 

194 in more detail. 

195 Comparison with RH relation in anisotropic plasmas: 

we We shall investigate the nature of these discontinuities formed along reconnection lay- 

19? ers by using RH relations in anisotropic plasmas [Karimabadi et al, 1995]. The basic 

us equations are as follows: 



WnW = 0, 



(10) 





pV n V t - e 



B n B t ] 1 



= 0, 



(12) 



4/T J 2 
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pV n V 2 + 



7 p\ . e + 2 T JB| 



- + 



-K. 



n i 



= 0, 



(13) 



J 2 



7 — 1 p J 3 47T 47T 47T 

where the brackets represent [A]\ = A ± — A 2 and the subscripts, 1 and 2, represent 
upstream and downstream, respectively. A n and A t respectively stand for the normal and 
tangential components of A. In above equations, the total pressure p = (p\\ + 2p±j /3 is 
introduced instead of using the double-adiabatic theory. (In more detail, see Karimabadi 
et al. [1995].) 

It is known that the modified intermediate Mach number M n = V n / (y/eVA n ) is useful 



to discuss the above RH solution, where Va u = B n /y/Airp is the Alfven velocity normal 
to the shock front [Hau and Sonnerup, 1989]. Then, the relation between the upstream 
modified intermediate Mach number M n \ and its downstream value M n2 are obtained, 
after some algebraic calculations, as follows: 



A a (e 2 , 6 U M A n2 ) ■ e\M A nl + 2A b (e 1 , e 2 , 6 U fa, M A n2 ) ■ e 1 M^ 1 + A c (e 1? e 2 , 6 u p u M' l n2 ) = 0, (14) 



where 



A„ = 



A b = 
A, = M, 



1 ^--^M^tan 2 ^, 

7 cos 2 0i 

"7-1 2(1- Cl ) , ei A 



7 



3 cos 2 2 cos 2 



- e 2 M 2 2 



„2 ) ' 2& 



7 + 1 



ei/?i 



7 



n2 



+ 



e 2 cos 2 ^i Ve 2 



+ ei£iM 2 2 tan 2 6 U 



' 2 7-2 _ 2 



and 

6 = 
6 = 



7-1 



M. 



7 

K 2 2 - 1) 



n2 
2 



7 



2 + i 

£2 



- tan 2 0! 



-e^tan 2 OA, 



1 / 1 

37 V e 2 
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211 0\ is the angle between the shock normal and the upstream magnetic field line. So four 

212 parameters, i.e., anisotropic parameters in both upstream and downstream regions, an 

213 upstream shock angle, and an upstream plasma beta are necessary to determine above 

214 RH relations. These parameters are calculated from physical quantities obtained by our 

215 simulation and we can obtain the relation between the upstream and downstream Alfven 

216 Mach numbers. Then, if one chooses a certain upstream Alfven Mach number, one obtains 

217 corresponding downstream Alfven Mach numbers. It should be noted that since plasma 

218 betas in two lobes are much smaller than unity, the upstream anisotropic parameter e\ is 

219 almost equivalent to unity as is shown in Figures 4-6. 

220 In Figure 7 we show such parameters in magnetic reconnection. From the top to the 

221 bottom, boundaries separating the upstream, transition, and downstream regions in the 

222 x-y plane, the inflow velocity, the upstream shock angle, and the downstream anisotropic 

223 parameter are shown, respectively. To eliminate the effects of initial current plasmas, 

224 we do not analyze the region x < —150 Ai. In Figure 7(a), we choose error bars for 

225 boundaries so that changes of both V x and V y are within them. Error bars for 9i, V n i, and 

226 e 2 stand for standard variations due to the averaging procedure for both upstream and 

227 downstream regions. In order to determine the upstream shock angle, we first evaluate 

228 the angle between the CPS, i.e., y — 0, and a shock surface, 9sc, by the method of the 

229 least square fit. Then, using magnetic field data obtained by our simulation results, we 

230 calculate the angle between the CPS and the upstream magnetic field. Finally, we obtain 

231 the upstream shock angle 9 1 , and the upstream velocity normal to the shock front, i.e., 

232 \V n i\ = \V x \sin9sc + \Vy\ cos @sc- As for the downstream anisotropic parameter, e 2 , we 



DRAFT 



January 23, 2012, 1:07am 



DRAFT 



X - 14 HIGASHIMORI AND HOSHINO: SLOW SHOCKS AND ION TEMPERATURE ANISOTROPY 

233 eliminated points whose standard variations are greater than 0.5, because of difficulties 

234 in identifying their shock downstream structures. 

235 Since Q\ varies over the range of Q\ ~ 78°-89° and the average of inflow velocity V n \ 

236 is nearly equal to 0.12 from Figures 7(b) and 7(c), we can suppose that the upstream 

237 Alfven Mach number normal to the shock surface, M nl = V n / (^/siVa cos 9 Br ^j , varies in 

238 the range M nl > 0.58 (M;^ > 0.33). As for the downstream anisotropic parameter, e 2 

239 increases with increasing the distance from the magnetic neutral point in Region 1, and 
2M e 2 is nearly equal to unity everywhere in Region 2. 

241 Now, we know e 2 , #i, and M nl . In addition, the upstream plasma betas are nearly equal 

242 to 10 -2 everywhere in the upstream region. Therefore, we can draw RH shock solution 

243 curves. Figure 8 shows such RH solutions in three cases: (/3i,6i,e 2 ) = (10~ 2 , 84°, 0.60), 

244 (10~ 2 , 84°, 0.85), and (10~ 2 , 84°, 1.00). Since there is no recognizable dependence of 6i on 

245 the nature of shock solution curves, Q\ is supposed to be 84° here. The curves in the 

246 region where M n \ < 1.0 correspond to slow shock solutions. As is shown in the earlier 

247 works [Lyu and Kan, 1986; Karimabadi et al, 1995], slow shock solutions are sensitive 

248 to the temperature anisotropy and the solution curves stretch as e 2 becomes large. From 

249 this figure, one can find that if e 2 is larger than 0.85, the minimum M 2 X value obtained by 

250 our simulation, i.e., M 2 X = 0.33 (M 2 : > 0.33), always has intersection points with the RH 

251 shock soolution curve and slow shock solutions can exist over the range 0.33 < M 2 X < 1.0. 

252 On the other hand in the case of e 2 < 0.85, the area, where slow shock solutions can exist, 

253 shrinks and it becomes harder for slow shocks to have their solutions with given Mach 
numbers. 
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255 Based on the above discussion, let us examine Figure 7(d) again. The horizontal dash- 

256 dotted line corresponds to e 2 = 0.85. In the region —115 A; < x, since e 2 is smaller than 

257 0.85, the range of upstream Alfven Mach numbers where slow shocks can exist becomes 

258 narrower. While in the outer region x < —115 A;, the ion temperature anisotropy is 

259 small and e 2 is nearly equal to unity. It enables slow shocks to exist stationarily in 

260 the calculated upstream Alfven Mach numbers. The compression ratio p 2 /pi in such 

261 regions [simulation/theory (error %)] can be calculated as [1.9/2.5 (32 %)] at x ~ —115 A; 

262 (Region 1) with M nl = 0.95, 6 X = 82.8°, ft = 10~ 2 , and e 2 = 0.98, and [2.1/2.6 (24 %)] 

263 at x ~ -145 A; (Region 2) with M nl = 1.0, 6 1 = 87.5°, ft = 10~ 2 , and e 2 = 0.95. 

4. Discussion and Conclusion 

264 We have discussed the relation between the ion temperature anisotropy obtained by 

265 our two-dimensional kinetic simulation of magnetic reconnection and the formation of 

266 slow shocks in accordance with the RH relations. From the point of view of RH relations 

267 in anisotropic plasmas, the parameters to determine shock solutions are the upstream 

268 plasma beta ft and shock angle 9 1 , the downstream ion temperature anisotropy e 2 , and 

269 the upstream Alfven Mach number M n \. Among these four parameters, ft, 6\, and 

270 e 2 have particularly an important influence on the nature of RH shock solution curves. 

271 We have evaluated the spatial profiles of these parameters in magnetic reconnection and 

272 discussed whether or not the discontinuities satisfy the conditions for slow shocks. In 

273 this study, it has been shown that the downstream ion temperature anisotropy along the 

274 reconnection layer decreases with increasing distance from the neutral point, and that a 

275 pair of current layers is formed in the region where plasmas are considerably isotropized 

276 {\x\ > 115 A;). This spatial distribution of the downstream ion temperature anisotropy is 
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277 strongly related to the formation of slow shocks, and the relaxation of the ion temperature 

278 anisotropy allows RH shock solutions in a broad range of upstream Alfven Mach numbers 

279 in collisionless magnetic reconnection. Let us discuss the dependence of such parameters 

280 on the formation of slow shocks in magnetic reconnection in more detail. 

281 First, we will refer to effects of downstream temperature anisotropy. Previous studies for 

282 the RH solutions in anisotropic plasmas suggest that the most important parameter for the 

283 RH relations is the downstream anisotropic parameter e 2 [Lyu and Kan, 1986; Karimabadi 

28 4 et al, 1995]. In collisionless magnetic reconnection, such temperature anisotropy is due to 

285 the PSBL beam ions, whose characteristic bulk velocity is about the lobe Alfven velocity. 

286 Plasma mixing between cold lobe plasma, whose bulk velocity is about a tenth of the lobe 

287 Alfven velocity, and PSBL ion beam components result in high temperature parallel to 

288 the magnetic field. Under these circumstances, the downstream anisotropic parameter e 2 

289 becomes smaller than unity. Here, we would emphasize that the generation mechanism of 

290 the ion temperature anisotropy in magnetic reconnection is quite different from that in the 

291 slow shocks without magnetic reconnection. In case of slow shocks without reconnection, 

292 incident ions and backstreaming ions from the shock downstream region are known to 

293 form the ion temperature anisotropy. On the other hand, in the case of slow shocks 

294 with magnetic reconnection, in addition to the backstreaming ions, PSBL ion beams 

295 accelerated around the diffusion region can contribute to the ion temperature anisotropy 

296 as well [Hoshino et al., 1998]. Getting back to the diagram of M^ 2 -M^ plot discussed 

297 before, as e 2 becomes small, the region where slow shock solutions exist shrinks. This is 

298 the main reason why slow shocks are hard to have their solutions near the diffusion region 

299 in magnetic reconnection. 
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00 In addition to the formation of such temperature anisotropy, a particular interest is why 

01 the relaxation of the ion temperature anisotropy occurs away from the diffusion region. In 

02 collisionless magnetic reconnection, the anisotropic PSBL ions can become thermalized by 

03 the Alfvenic waves generated by the ion-cyclotron beam instability, the ion/ion cyclotron 

04 instability [Winske and Leroy, 1984], and the EMIIC (electromagnetic ion/ion cyclotron) 

05 instability [Winske and Omidi, 1990]. From observations in the Earth's magnetotail, the 
oe existence of the ion cyclotron beam instability in the PSBL has been confirmed [Kawano 

07 et al, 1994; Takada et al, 2005]. Additionally, the fire-hose instability can take place in 

08 the reconnection exhaust [Karimabadi et al, 1999; Liu et al, 2011]. Our simulation results 

09 also support the fire-hose instability enhanced near the CPS. As shown in Figures 4-6, 

10 since the temperature anisotropic parameter e is greater than 0, the fire-hose instability 
a is absent in the transition region. However, near the CPS, the condition e < for the 
12 fire-hose instability can be satisfied in —100 A ; < x. We think that these instabilities 
io may play an important role in the relaxation process of the downstream ion temperature 
14 anisotropy. As a result, Maxwellian-like distribution functions of ions are observed away 
is from the neutral point. From the viewpoint of anisotropic RH relations, such relaxation 
i6 of the downstream temperature anisotropy away from the diffusion region enables slow 
i? shocks to have their solutions in wider range of the upstream Alfven Mach numbers. 

is We briefly refer to the upstream shock angle, 6\. From our simulation results, Q\ is 

19 more or less in the range of 78°-89°. According to the anisotropic RH relations, such a 

20 variation obtained by our simulation does not result in a recognizable impact on the nature 

21 of solution curves. However, if 9\ is smaller than 78° by keeping the inflow velocity of 0.12, 
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322 the upstream Alfven Mach number becomes in the range of M n \ < 0.58 (M^ 1 < 0.33) and 

323 the slow shock condition is not satisfied. 

324 In this paper, we have shown the existence of slow shock discontinuities and the resul- 

325 tant current sheet profiles in both Regions 1 and 2, but the behavior of the downstream 

326 temperature anisotropies are different between two regions. The downstream ion temper- 

327 ature in Region 2 is nearly isotropic, while that in Region 1 is rather high even though the 

328 temperature anisotropy can be relaxed with increasing distance from the X-type neutral 

329 point. This behavior seems to be important for the satellite observation of the slow mode 

330 shock. So far the observational study of the slow shock detection in the Earth's magneto- 

331 tail assumed the isotropic temperature [Saito et ai, 1995; Seon et ai, 1995], but the RH 

332 study of the slow mode shock including the temperature anisotropy may distinguish the 

333 slow shock region between Regions 1 and 2. 

334 In Region 1, we found the relaxation of the anisotropic temperature and the formation 

335 of the slow shock, but one might indicate that larger scale simulations would result in 

336 more isotropic plasma distribution in the slow shock region. In fact, we studied larger 

337 scale simulations, but we found that as time goes on, other magnetic islands are formed 

338 from the diffusion region, grow, and are ejected into the outflow region. As a result, the 

339 size of Region 1 cannot become larger than ~ 120 Ai. 

340 We obtained that the minimum distance for relaxation of the temperature anisotropy 

341 required for the formation of slow shock is about 115 times ion inertia length based on our 

342 two-dimensional hybrid simulation. However, some other processes that are not included 

343 in our simulation may quickly reduce the ion temperature anisotropy in a shorter spatial 

344 scale from the X-type neutral point. Such candidates might be the three dimensionality 
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3« (e.g., drift-mode instabilities), or the instabilities due to electron kinetic effects which 

346 cannot be considered in our hybrid model. Recently Yin et al. [2007] performed the 

347 oblique slow shocks by full-particle simulations, and discussed how electron kinetic effects 

348 are related to kinetic Alfven waves and could alter the structure of slow shocks. These 

349 possibilities will be investigated in future works. 
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Figure 1. (a) Magnetic field lines and (b) the electric current in the z direction at time 
t = 415 Of . The initial current is in the — z direction. 
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Figure 2. Enlarged views of the reconnection layer at x < 0. (a) Magnetic field lines, 
(b) the out of plane magnetic field B z , (c) the ion flow vector, (d) the ion temperature, 
and (e) the ion temperature ratio Ty/T^i. are shown in the x-y plane. 
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Figure 3. Ion velocity distribution functions (f(v x ,v y )) in the upstream, transition, 
and downstream regions. Their locations are indicated as white squares in Figure 2(e). 




Figure 4. Cross-sectional views of discontinuities at x = —88.3 Ai. 
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Figure 5. Cross-sectional views of discontinuities at x — —115.0 A ; . 




Figure 6. Cross-sectional views of discontinuities at x — —145.0 A ; . 
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Figure 7. (a) Dash-dotted line and dotted line respectively stand for the boundaries 
separating the upstream and transient regions, and the ones separating the transient and 
downstream regions. In addition, spatial profiles of (b) the upstream shock angle, (c) 
the upstream velocity normal to the surfaces of discontinuities, and (d) the downstream 
anisotropic parameter are shown as a function of the distance from the neutral point. 
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Figure 8. M^ a -M^ lX plot of RH solutions. Isotropic (e 2 = 1) and anisotropic cases 
(e 2 = 0.85, 0.6) are shown. The upstream plasma beta and the shock angle are assumed 
to be (3i = 10~ 2 and 6i = 84° in all these three cases. The square of the minimum 
upstream Alfven Mach number calculated in our simulation, i.e., = 0.33, is shown 
by the horizontal dotted line. The area where slow shocks (SS) can exist is indicated by 
a double-headed arrow (0.33 < M%i < 1.0 in case of e 2 > 0.85). 
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